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Abstract 

We revisit the short-time dynamics of 2D Ising model with three spin 
interactions in one direction and estimate the critical exponents z, 0, (3 
and v. Taking properly into account the symmetry of the Hamiltonian 
we obtain results completely different from those obtained by Wang et 
al.. For the dynamic exponent z our result coincides with that of the 
4-state Potts model in two dimensions. In addition, results for the static 
exponents v and /3 agree with previous estimates obtained from finite size 
scaling combined with conformal invariance. Finally, for the new dynamic 
exponent 6 we find a negative and close to zero value, a result also expected 
for the 4-state Potts model according to Okano et al.. 



Since the work by Janssen et al ^ and Huse [|| pointing out the existence 
of another universal stage at an early time critical dynamic, several statistical 
models have been investigated to confirm the analytical predictions about the 
"critical initial slip" and to enlarge the knowledge of critical phenomena 
[il' H' @' 0' @- "^^^ investigation of the universal behavior in short-time 
dynamics avoids the critical slowing down effects of the equilibrium and provides 
an alternate way |^ to calculate the new exponent 9 which governs the behavior 
of the magnetization, the dynamic critical exponent z as well as the static 
exponents f3 and i^. 

In this letter we adopt this approach to study the two-dimensional (2D) 
Ising model with three-spin interactions in one direction and calculate its set of 
exponents. Motivation came from the fact that in a recent paper by Wang et al 
[Eol estimates obtained were in complete disagreement with pertinent results. 



Here we show that, when the symmetry of the model is taken properly into 
account, good agreement is obtained with expected results. 

The Hamiltonian of the 2D Ising model with three spin interactions (m = 3) 
in one direction is pl| 

{m-l 
KyS,,,S^+i.j+K,Y[S,,j+A (1) 
1=0 J 

where Si^ = ±1 is the Ising spin variable. The model is known to be self-dual 
[ p!2| , its critical line being 



sinh 2Kr, sinh 2Ky = 1 (2) 
for all m. For the particular isotropic case (K^ = Ky) the critical coupling is 

= J/keTc = i ln(V2 + 1) = 0.4406867, (3) 

which is the same of the standard 2D Ising model. 

Symmetry analysis and ground state degeneracy considerations suggest that 
the model is in the same universality class as the q-state Potts model, whenever 
q = 2™^^ Q. This result is supported by finite-size scahng studies [14|, 
[12 , pq| , weak and strong coupling expansio ns pq |, conformal invariance |17|, 
[18 , ]19|, standard Monte Carlo simulations [pO[, [pl[[[p^ and mapping of the 
TO = 3 model in the extreme anisotropic limit of the 4-state Potts model [^ . In 
most of those papers the argument to include the Ising model with three-spin 
interaction and the 4-state Potts model in the same universality class is based 
on the value of the exponents v and q;(w 2/3) [Q. All of them respect the 
symmetry of the Hamiltonian. Very little is known about the exponents (3 [^ 
and z Q. 

The ground state for the general TO-spin interaction is 2™~^degenerated [Q, 
which implies that the 2D Ising model with to = 3 spin interaction is four- 
fold degenerated. The relevant symmetry of this model is semi-global [|2^ and 
the Hamiltonian is symmetric under the reversal of all the spins in any two 
sublattices, which leads to the existence of three independent interpenetrating 
sublattices in the system. At T = 0, the possible states consist of repetitions of 

the patterns H , — I — , h in the horizontal direction, copied along 

the lines. According to preceding statement, it is important to take lattice sizes 
that are multiple of three - in order to respect the symmetry of the Hamiltonian - 
and to work with the appropriate order parameter - the sublattice magnetization 
- in order to avoid the effect of staggered magnetization. We suspected that 
Wang et al [ pO[ didn't consider the m-spin symmetry in their simulations neither 
worked with the magnetization of sublattice, since they just presented results 
for square lattices which do not obey that previous condition (multiple of three). 
In this sense, it seemed to us relevant to repeat simulations, paying attention 
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to the above mentioned points, to check the apparent failure of the short-time 
approach in this case. 

We began by repeating the analysis made by Wang et al Q for the Binder 
cummulant 

U^UL)JJf^-, (4) 
(M(f)) 

L L 

where < > means average on samples, t is the time, M (t) = J2 J2 ^ij the 

i=ij=i 

magnetization at time t and L the size of the square lattice. They argue that 
this expression obeys the power law form: 

C/(t, L) a (5) 

when the dynamical process starts from an ordered state (mo = 1), which is a 
fixed point under renormalization group transformation. Fig. 1 shows explicitly 
the different results obtained when we use the magnetization of the sublattice. 
The average is taken over 50000 independent initial configurations and the error 
bars (smaller than the size of the points) are obtained by repeating five times 
each simulation. When the simulation is performed without the sublattice con- 
siderations the slope of the curve agrees with that presented by Wang et al and 
confirmed our suspicions. From our point of view | |2^ , though, this cummu- 
lant should obey the power law form (|^) only when different initial conditions 
were used in the study of the magnetization and its second moment. Scaling 
arguments |^ assert that the second moment of magnetization behaves as 

(Af2(t)) 0C<('^-2/5/-)/-, (6) 

only when the samples are taken with zero initial magnetization (mg = 0). On 
the other hand, short-time scaling behavior implies 

{M{t)) a t-^/"'' (7) 

for samples starting from the ordered state (mg — 1) Thus, performing 

two different simulations under those mentioned conditions, we obtain the time 
evolution of the ratio (A'/^(i)) / (Af(t))^ which furnishes the exponent d/z in a 
log-log plot (Fig. 2). From the slope of that curve we estimate z = 2.380 ± 
0.004. 

To confirm our result we used two other approaches: the generalized fourth 
order Binder's cumulant |^ and the parameters Q and R introduced by de 
01iveira|p7|. In both cases, the collapse of the curves for different lattice sizes, 
at critical temperature, are used to determine the dynamical critical exponent 
z from short-time simulations. The Binder's cummulant, 

Ui{t,L) = l--^ ^ (8) 

3(M(2))' 
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satisfies, at T = Tc the scaling relation 



Ui{t,L) = U4ib'H,b-^L) (9) 

where b — L/L' since U4 scales as L°. This technique has proved to be useful 
in determining the exponent z and was applied to the 2D and 3D Ising models 
[||, |28j, the 3-state Potts model 0, the majority vote model |^ and cellular 
automata The initial magnetization of samples, in this case, is zero as well 
as the correlation length. The results are good enough. Error bars, however, 
are bigger than those obtained by damage spreading technique [Q. In Fig. 3, 
we show the collapse of the cumulant U4 for lattice pairs {L, 2L) with z = 2.3. 
In fact, the range of z for which the collapse is still observed is given by 
z — 2.3 ± 0.1. This result supports our previous estimate for z (2.383 ± 0.004) 
and can be related to the 4-state Potts model exponent [|2| . In order to stress 
the importance of considering the symmetry of Hamiltonian, we exhibit in Fig. 
4 the deformation of the Binder cumulant when there is no sharp preparation 
of the initial magnetization on the sublattices and the magnetization evolution 
is calculated without restrictions. 

When we use the parameters Q{t,L) — l^sign {^j2M{t)')'^ and R{t,L) = 

(( sigrij^ J2 ^i] \ sign J- J2 ] ) °f '^^ 01iveira]27t and scaling relations for 
V bottom J \ top j 

T^TrM 



Qit,L)^Qib-H,b-^L) 

and 



R{t,L) ^ R{b-H,b-^L) 

we obtain the collapse among curves (see Fig. 5 and Fig. 6) for different lattices 
when time is scaled with z — 2.3. Samples in these cases were initialized with 
all spins up (toq = 1). 

In order to calculate the exponent we studied the derivative of the mag- 
netization, which presents the following scaling form: 



dr\nM{t,T)\r=0^t^'''''dr'\nF{T')\,,=0. (10) 

Fig. 7 shows the power law behavior of the 9i- In M(f,r) when Ar = 0.0002. 
The measured slope of the curve gives 1/vz — 0.624 ± 0.005. Thus, taking 
z = 2.383 ±0.004 we find v = 0.67 ±0.01, to be compared with the conjectured 
value V = 2/3. 
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Since the values of the exponents v and z are already known, the exponent /3 
of the magnetization can be obtained from the power law increase of the second 
moment of magnetization, Eq.(||). Fig. 8 presents, in a double- log scale, the 
polynomial behavior of (^M'^'^\t)). From the slope of these lines we estimate 
p = 0.11 ± 0.02. This result is in agreement with the expected value /3 = 0.125 

0, m. 

In order to extend the picture of universality we investigate the exponent 9 
by a recent technique proposed by Tome and de Oliveira |36|. In their paper 
they show that the exponent 9 can also be independently calculated in despite 
of the sharp preparation of the samples, since the time correlation of the to- 
tal magnetization in samples with a random initial configuration also exhibits 
polynomial behavior 

(Af(t)Af(O)) = 1 ^EE^«(*)^^(0)) « (11) 

This procedure avoids the use of an initial state with a nonzero (but small) 
magnetization as well as the numerical extrapolation mo — > to calculate the 
dynamic exponent 6. Fig. 9 shows, in a log-log scale, the (Af (i)Af(O)) behavior 
for different lattice sizes. The slope of those curves give us ^ = —0.03 ± 0.01 
which is compatible with the conjecture by Okano et al for the 4-state Potts 
model §. 

In summary, we have obtained static and dynamic critical exponents for 
the Ising model with multispin interactions using short-time Monte Carlo sim- 
ulations. Our results show that this model and the 4-state Potts one share 
the same set of critical exponents even at dynamic level. When compared to 
the paper by Whang et al this letter shows the importance of taking properly 
into account the symmetry of the Hamiltonian to deal with magnetization and 
boundary conditions. We stress that the present result for the exponent v is 
better than previous estimates obtained by finite size scaling and Monte Carlo 
approach Q. 
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Figure 1 : Log-log plot of the cumulant U=<M^>/<M>^-1 versus time. Our results were 
obtained using the sublattice magnetization (upper curve), but samples (50000) were 
originally in the ordered state, as done in the paper by Wang et al [10] (lower curve). 
The error bars are smaller than the symbols. 




Figure 2: Log-log plot of the ratio L/=</W^>/</W>^ when different initial conditions 
are used to measure the magnetization {m^=^) and its second moment (mQ=0). 
Error bars were calculated over 5 sets of 50000 samples and are smaller than 
the dot sizes. 
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Figure 3: Generalized Binder cumulant for L=9, 18,36 (continuous curves) compared to 

the rescaled cumulant in time through t-t/2' for L=18,36, 72 (dotted curves). The best fit 
is obtained when z=2.3. 




time 



Figure 4: Deformation of the generalized Binder cumulant {L=12) when the sublattices 
are not considered in the Ising model with three spin interaction. 
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Figure 5: Time behaviour of Q(t) for lattices L=18, 36, 48 (full lines) and the 
corresponding time rescaled curves for lattices L=12, 24, 36 (open circles). 
The collapse was obtained with z=2.3. At T=0 all the spins are in the up 
direction {nig=1). Statistical errors are smaller than the symbols. 
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Figure 6: Short-time behaviour of R(t) of lattices L=18, 24, 36 (full lines) together 

with the corresponding time rescaled curves for the lattices L=24, 36, 48 (open symbols). The 

initial condition is given by m^=1 and error bars are smaller than the symbol 

sizes. The best fit is obtained with z=2.3. 




Figure 7: Power law behaviour of the derivative of the magnetization 
(eq. 12) gives the critical exponent v of the correlation lenght. 




Figure 8: Log-log plot of the temporal increase of the second moment of the 
magnetization. The error bars were calculated over 5 sets of 50000 samples 
and are smaller than the dotsizes. The initial magnetization is given by m =0 




Figure 9: Time correlation of the total magnetization for samples which have 
in average 111^=0. Error bars were calculated over 5 sets of 50000 samples. 



